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Abstract. We determine the universal law for fidelity decay in quantum computations of complex dynamics 
in presence of internal static imperfections in a quantum computer. Our approach is based on random 
matrix theory applied to quantum computations in presence of imperfections. The theoretical predictions 
are tested and confirmed in extensive numerical simulations of a quantum algorithm for quantum chaos 
in the dynamical tent map with up to 18 qubits. The theory developed determines the time scales for 
reliable quantum computations in absence of the quantum error correction codes. These time scales are 
related to the Heisenberg time, the Thouless time, and the decay time given by Fermi's golden rule which 
are well known in the context of mesoscopic systems. The comparison is presented for static imperfection 
effects and random errors in quantum gates. A new convenient method for the quantum computation of 
the coarse-grained Wigner function is also proposed. 



PACS. 03.67.Lx Quantum Computation 
Quantum chaos; semiclassical methods 



1 Introduction 
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a search of unstructured database which has a quadratic 
speedup comparing to any classical algorithm 



Recently a great deal of attention has been attracted to 
the problem of quantum computation (see e.g. [llllllS|)» A 
quantum computer is viewed as a system of qubits. Each 
qubit can be considered as a two-level quantum system, 
e.g. one-half spin in a magnetic field. For Uq qubits the 
whole system is characterized by a finite - dimensional 
Hilbert space with N = 2^*? quantum states. It has 
been shown that all unitary operations in this space can 
be realized with elementary quantum gates which include 
one-qubit rotations B'^^^ and two-qubit controlled oper- 
ations, e.g. controlled-NOT gate C^^^ or controlled 
phase-shift gates B^'^^cj)) (see e.g. 3,4 ). The gates C^^^ 
and B^'^\(j)) assume that the interaction between qubits 
can be switched on and off in a controllable way with suf- 
ficiently high accuracy. Various computational algorithms 
in the space N can be represented as a sequence of ele- 
mentary gates. A general unitary operation (unitary ma- 
trix) in this space requires an exponential ( in n^) number 
of elementary gates. However, there are important exam- 
ples of algorithms for which the quantum computation can 
be performed with a number of operations (gates) much 
smaller than with the classical algorithms. The most fa- 
mous is the Shor algorithm for factorization of integers 
with Uq digits which on a quantum computer can be per- 
formed with O(n^) gates contrary to an exponential num- 
ber of operations required for any known classical algo- 
rithm 5 . Another example is the Grover algorithm for 



A quantum computation can be much faster than a 
classical one due the massive parallelism of many-body 
quantum mechanics since any step of a quantum evolu- 
tion is a multiplication of a vector by a unitary matrix. 
A very important example is the quantum Fourier trans- 
form (QFT) which can be performed for a vector of size 
N = 2"^^ with 0{nq^) gates instead of 0(ng2^«) classical 
operations required for the fast Fourier transform (FFT) 
(see e.g. PISl)- With the help of QFT the quantum evolu- 
tion of certain many-body quantum systems can be per- 
formed in a polynomial number of gates [7||H|. Another 
example can be found in the evolution of quantum dy- 
namical systems which are chaotic in the classical limit 
(see e.g. |9jil0|). Such systems are described by chaotic 
quantum maps and include the quantum baker map JT], 
the quantum kicked rotator the quantum saw-tooth 
map 13 and the quantum double- well map 14 . For them 
a map iteration can be performed for A^-size vector in 
Oijiq^) or Oijiq^) gates while a classical algorithm would 
need 0(ng2^'?) operations. This however does not neces- 
sary lead to an exponential gain since the final step with 
extraction of information by measurements also should be 
taken into account. Thus, for example, the quantum sim- 
ulation of the Anderson metal-insulator transition gives 
only a quadratic speedup even if each step of quantum 
evolution is performed in a polynomial number of gates 
P^j . Among other algorithms, let us refer to the quantum 
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computation of classical chaotic dynamics where some new 
information can be obtained efficiently p"6nT7) . 

The main obstacle to experimental implementation of 
a quantum computer is believed to be decoherence induced 
by unavoidable couplings to external world (see e.g. TS^). 
However, even if we imagine that there are no external 
couplings there still remains internal static imperfections 
inside a quantum computer. These static imperfections 
generate residual couplings between qubits and variation 
of energy level-spacing from one qubit to another. As it 
was shown in 19 such imperfections lead to emergence 
of many-body quantum chaos in a quantum computer 
hardware if a coupling strength exceeds a quantum chaos 
threshold. In a realistic quantum computer this thresh- 
old drops only inversely proportionally to the number of 
qubits Uq while the energy spacing between nearby lev- 
els drops exponentially with Uq. The dependence of this 
threshold on quantum computer parameters was studied 
analytically and numerically by different groups 19 ,20 ,21 , 
E^lESj The time scales for onset of quantum chaos were 
also determined. 

It is of primary importance to understand how effects 
of external decoherence and internal static imperfections 
affect the accuracy of quantum computations. A very con- 
venient characteristic which allows to analyze these ef- 
fects is the fidelity / of quantum computation. It is de- 
fined as f{t) = I < tlj^{t)\tlj{t) > p where > is 
the quantum state at time t computed with perfect (or 
ideal) gates, while \^l^e{t) > is the quantum state at time 
t computed with imperfect gates characterized by an im- 
perfection strength e. If the fidelity is close to unity then 
a quantum computation with imperfections is close to the 
ideal one while if / is significantly smaller than 1 then the 
computation gives, generally, wrong results. 

At first the fidelity was used to characterize the ef- 
fects of perturbation on quantum evolution in the regime 
of quantum chaos j23j- Indeed, for the classical chaotic 
dynamics the small errors grow exponentially with time 
while for the quantum evolution in the regime of quantum 
chaos small quantum errors only weakly affect the dynam- 
ics. For example, the time reversibility is broken by small 
errors for classical chaotic dynamics while it is preserved 
for the corresponding quantum dynamics in presence of 
small quantum errors 25 ,26 . In the context of quantum 
computation the qualitative difference between classical 
and quantum errors is analyzed in ^^l- Recently, the in- 
terest to the fidelity decay induced by perturbations of 
dynamics in the regime of quantum chaos has been re- 
newed 27,28,29,30,31,32,33,34 . It has been shown that 
the rate F of exponential decrease of / is given by the 
Fermi golden rule for small perturbations while for suffi- 
ciently strong perturbations the decay rate is determined 
by the Kolmogorov- Sinai entropy related to the Lyapunov 
exponent of classical chaotic dynamics 2 8. 29^ . For small 
perturbations the fidelity decay can be expressed with the 
help of correlation function of quantum dynamics that al- 
lows to understand various peculiarities of the decay 31 . 

Until recently the fidelity decay and accuracy of quan- 
tum computations have been mainly analyzed for the case 



of random noise errors in the quantum gates 135111611131 
[1^136113 7j . Quite naturally in this case the rate of fidelity 
decay is proportional to the square of error amplitude s 
{r (X e^). Indeed, a random error of amplitude e trans- 
fers a probability of order from the ideal state to all 
other states and as a result the fidelity remains close to 
unity (within e.g. 10% accuracy) during a time scale tf ^ 
l/{e^ng). Here rig is the number of gates per one map it- 
eration and for polynomial algorithms Ug ~ {e.g. for 
the quantum saw-tooth map 7 = 2 13 ). 

Contrary to the case of random errors the effects of 
static imperfections on fidelity decay have been studied 
only in 13736' . The numerical simulations performed there 
with up to 18 qubits show that for small strength of static 
imperfections e the time scale varies as ~ l/{eng.^/n^). 
Such a dependence implies that in the limit of small s the 
effects of static imperfections dominate the fidelity decay 
comparing to the case of random errors 13,36 . Simple es- 
timates based on the Rabi oscillations have been proposed 
to explain the above dependence extracted from numerical 

data nnusni- 

Since the numerical results show that the static imper- 
fections lead to a more rapid fidelity decay, compared to 
random errors fluctuating from gate to gate, it is impor- 
tant to investigate their effects in more detail. This is the 
aim of this paper in which we carry out extensive numeri- 
cal and analytical studies of static imperfections effects on 
fidelity decay using as an example a quantum algorithm 
for the quantum tent map which describes dynamics in 
a mixed phase space with chaotic and integrable motion. 
For the case when the algorithm describes the dynam- 
ics in the regime of quantum chaos a scaling theory for 
universal fidelity decay is developed on the basis of the 
random matrix theory (RMT) ^38|39,4Q^. This theory is 
tested in extensive numerical simulations with up to 18 
qubits and the obtained results confirm its analytical pre- 
dictions which are rather different from the conclusions of 
the previous studies 13,36 . We also investigate the regime 
of fidelity decay for integrable quantum dynamics where 
the situation happens to be more complicated. In addi- 
tion, a simple quantum algorithm is proposed for approx- 
imate computation of the coarse-grained Wigner function 
(the Husimi function) j^TlE2| and its stability in respect 
to imperfections is tested on the example of quantum tent 
map. 

It is important to note that all quantum operations re- 
quired for implementation of the quantum tent map have 
been already realized for 3 - 7 qubits in the NMR-based 
quantum computations reported in Refs. 43,44 . An effi- 
cient measurement procedure for fidelity decay in quan- 
tum computations is proposed in 45 . 

The paper is organized as follows. In Sec. 2 we describe 
the classical and quantum tent map. The algorithm for 
quantum dynamics is derived in Sec. 3. The fidelity decay 
for random errors in quantum gates is analyzed in Sec. 4. 
The analytical theory for fidelity decay induced by static 
imperfections is developed on the basis of RMT approach 
in Sec. 5. This theory is tested in extensive numerical sim- 
ulations presented in Sec. 6. An approximate algorithm for 
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the quantum computation of the Husimi function is stud- 
ied in Sec. 7. The conclusion is given in Sec. 8. 



2 Classical and quantum tent map 

We consider a kicked rotator whose dynamics is governed 
by the time dependent Hamiltonian, 
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with the potential of kick 



v{e) 



{0 - n){0 - 



27r) 



5(1 - m) 



0<0 <n 
TT < 6* < 27r 



(1) 



(2) 



where is taken modulo 27r and 5{t) is a (^-function, m is 
an integer. The parameter k determines the kick strength 
and T gives the rotation of phase between kicks. It is easy 
to see that the classical evolution for a finite time step 
t ^ t -\- 1 with respect to the Hamiltonian (QJ can be 
described by the map, 



p = p 



V\0) 



= 0^pT { mod 27r) 



(3) 



Here bars mark new values of the dynamical variables after 
one map iteration. This map is similar in structure to the 
Chirikov standard map |42|. The derivative of the kick- 
potential, 



0) 



9) 



< < TT 

TT < 6* < 27r 



(4) 



has a tent form and is continuous but not differentiable at 
^ = and = 71. This is an intermediate case between the 
standard map 0^1 with a perfectly smooth kick-potential 
and the saw-tooth map 13* with a non-continuous poten- 
tial. 



k = 0.53 



fc = 4/3 



k = 1.7 




Fig. 1. Classical Poincare sections of the map © in (6*, p) 
plane for T = 1 and i^^ = /c = 0.53, K = k = 4/3 and K = k = 
1.7. 



K the dynamics is governed by a KAM-scenario with the 
Kolmogorov-Arnold-Moser (KAM) invariant curves and a 
stable island at^ = 37r/2,p = and a chaotic layer around 
separatrix starting from the unstable fixed point (saddle) 
at ^ = 7r/2, p = 0. At K = 4/3, the last invariant curve 
is destroyed and one observes a transition to global chaos 
with a mixed phase space containing big regions with reg- 
ular dynamics |T7lllH) . 

In Fig. n the Poincare sections of the map © for the 
three values K = 0.53, 4/3, 1.7 are shown. Here we have 
replaced p by its value modulo 27r/T which is appropri- 
ate since the classical map is invariant with respect to the 
shift p ^ p -\- 27r/T. Fig. confirms the above scenario of 
a transition to global chaos at K = 4/3. In the following, 
we are particularly interested in a typical case K = 1.7, 
which exhibits global chaos with quite large stable islands 
in phase space related to the main and secondary reso- 
nances. 

At > 4, the phase space becomes completely chaotic 
and the dynamics is characterized by a diffusive growth in 
p. The diffusion rate D inp can be obtained with the help 
of random phase approximation that gives 



{{p-po)')^Dt , D = {V' {9)^)0 



12 



(5) 



Here and below t is an integer which gives the number of 
map iterations (kicks). 

The quantum dynamics of the Hamiltonian Q is given 
by the Schrodinger equation, 

i^^\m> = Hit)\m> ■ (6) 

Here in the Hamiltonian H{t) the variables p and are 
operators with the commutator [p,0] = —i. They have 
integer eigenvalues p for p and real eigenvalues in the 
interval [0,27r[ for 0. As in the classical case one can de- 
termine the evolution for one map iteration: 

\^{t + 1)> = /7 \^{t)> = e-'^p"/^ e-'^^^^ \^{t)> . (7) 

Eq. {TJ corresponds to the quantized version of the clas- 
sical map © and defines a quantum map that can be 
efficiently simulated on a quantum computer. Here h = 1 
and the quasiclassical limit correspond to T ^ 0, /c ^ oo 
with K = kT = const. 

For quantum dynamics we concentrate our studies on 
the case K = kT = 1.7 and T = 2i\ jN that corresponds 
to evolution on one classical cell (see Fig.^ with N quan- 
tum states. As an initial state we use a minimal coherent 
wave packet corresponding to a given N which is placed 
in a chaotic or integrable component (near the unstable 
or stable fixed point at p = 0,^ = 7r/2 or ^ = 37r/2). An 
example of the Husimi function (see Sec. 7) for a chaotic 
case is shown in Fig. |2|for different moments. 



The dynamics of the classical tent map (|2| depends 
only on one dimensionless parameter K = /cT, its prop- 
erties have been studied in j47ll48j . For small values of 



3 Quantum algorithm 

The quantum map Q can be simulated on a quantum 
computer in a polynomial number of gates. The quantum 
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Fig. 2. Density plots of the Husimi functions (see below Sec. 7) 
of the quantum state \\\){f)> for t = 5 (top left panel), t = 15 
(top right panel) and t — 5625 (bottom left panel) with an 
initial state |'0(O)> chosen as a minimal coherent (gaussian) 
wave packet closely located to the saddle at ^ = 7r/2, p = 0. 
The bottom right panel corresponds to the quantum state after 
5625 iterations computed on a quantum computer with static 
imperfections for £ = 7 • 10~^, here the fidelity is / = 0.9388 
(see below Sees. 3,6 for details). The density is minimal for 
blue/black and maximal for red/white. All panels correspond 
to Uq — 16 qubits, i.e. to a finite dimensional Hilbert space of 
dimension AT = 2^^, T = ^i^jN. 

algorithm has similarities with those described in p"2HT3| . 
To perform one iteration of the quantum map Q we repre- 
sent the states |?/;> by a quantum register with qubits. 
In particular, the eigenstates \p> of the momentum oper- 
ator are identified with the quantum-register states 

|aO, Q^l, . . . , Oin^-\> = \ao>0 |ai>l . . . |Q^ng-l>ng-l (8) 

where aj = or 1 and 

riq — l 

p=^a,2^. (9) 

j=o 

The states \0>j and |l>j correspond to the two basis 
states of the j— th qubit. Obviously, this representation 
introduces a Hilbert space of finite dimension N = 2^«; 
the operator p has the eigenvalues: p = 0,...,7V — 1. 

A quantum computer is a machine that is able to pre- 
pare a quantum register with a well defined initial condi- 
tion and to perform certain well controlled unitary oper- 
ations on this quantum register. These particular opera- 
tions are called quantum gates and one typically assumes 
that the quantum computer can be constructed with quan- 
tum gates that manipulate at most two qubits. Here we 



use as elementary gates the phase-shift gates and 
controlled phase-shift gates B^jJ : 

B^^\^)\p> = e'''^^\p> , (10) 

^S^W|P> = e^"^"'^|P> , jy^k (11) 

where p is of the form Q. These gates provide a phase 
factor e'^^ if aj = 1 for the simple phase-shift or if aj = 
a/e = 1 for the controlled phase-shift. Using Eq. Q, one 
easily verifies that the momentum dependent factor of the 
unitary operator in Q can be expressed in terms of 
these gates by: 

riq — l Hq — 1 

g-iTpV2^ -Q Bf^{-T2^+'') W Bf\-T2'i-^) . 

j<k j=0 

(12) 

The situation is different for the phase factor containing 
the kick-potential since this factor is not diagonal in mo- 
mentum representation. It is therefore necessary to trans- 
form to the basis of eigenstates of the operator 0. For this, 
following 1 we consider the unitary operator Uqyt de- 
fined by: 

AT-l 

[/QPT|p> = -=^e2-^^~/^|p>. (13) 

p—O 

Then the eigenstates of with eigenvalues = ^jf- are 
simply given by: t/qp^ \p> and more generally the opera- 
tors 6 and p are related by: 

^ = t^QFT (^) f^QFT. (14) 

Using Eqs. ((2j, Q and (|14|) it is straight-forward to show 
that the unitary factor of U containing the kick-potential 
can be written as: 

nq-2 

j<k 

X Bg)(fc7r2 2^+'=-2"^+2)} (15) 

nq-2 

X n {C-i(-^^'2^--^+2(2^--^+^ - 1)) 

X B^^\k7r^2^-''^{2^-''^^^ - 1))} Uqft • 

Here we have used a three-qubit gate for a controlled- 
controlled phase-shift defined by (cf. Eqs. (PT)|) . ([TT|) ): 

BjkiW \P> = e'"^-"^"^^ \p> . (16) 

In Eq. ([T5|) it appears for / = — 1 because of the two dis- 
tinct cases in Eq.©. In principle, this gate is not directly 
available in the set of one- and two-qubit gates used for a 



K.M.Frahm, R.Fleckinger and D.L.Shepelyansky: Quantum chaos and random matrix theory for fidelity decay 



5 



quantum computer. However, it can be constructed from 
5 two-qubit gates by: 



B 



(3)/ 



B 



(2) 



- 1 B^^^ 
2 / 



where C|.[^^ is the controlled-not gate that exchanges the 
states |0 >fe and |1 >k for the k-th qubit if o;; = 1. In 
matrix representation it is given by: 



(17) 



r 



(AT) 
kl 



/I 0^ 
10 
1 

\0 1 0, 



(18) 



where the index ai corresponds to the outer block struc- 
ture and ak to the inner block structure. 

Following the description 1 , the QFT operator [/qft 
can be written in the form: 



(19) 



where R is the unitary operator that reverses the order 
of the qubits and Aj is a one-qubit gate with the matrix 
representation: 



1 1 
1 -1 



(20) 



We note that in the outer product in Eq. (|T^ the factors 
are ordered from left to right with increasing j. 

Combining, Eqs. ([T^ . (P3|) and (|19|) . we see that the 
quantum map {TJ can be expressed by a total number of 
rig = — + 4 elementary quantum gates (and 2 R- 
operations). On a classical computer one iteration of the 
quantum tent map requires 0(ng2^'?) operations coming 
mainly from the FFT. 

To investigate the stability of the quantum algorithm 
for the tent map we consider two models of imperfections. 
The first model represents the random errors in quantum 
gates fluctuating in time from one gate to another (ran- 
dom noise errors). In this case for all phase-shift gates 
we replace (j) hy <p -\- Scj) with random S<p G [— £, e] that 
is different for each application of the gate. For the gates 
containing the Pauli matrix ax we replace it by n-cr where 
n is a random unit- vector close to Cx with \n — ex\ < s. 

The second model describes only static imperfections 
and is similar to one used in 19,13,36,15 . In this case the 
effect of static imperfections is modeled by an additional 
unitary rotation between two arbitrary gates which has 



the form : Us 
by: 



. Here the Hamiltonian SH is given 



(x) (x) 



a, 



(21) 



j=0 



where <jj^^ are the Pauli matrices acting on the jth qubit 
and Sj^ Jj G [— V^e, V^e] are random coefficients which 



are drawn only once at the beginning and kept fixed during 
the simulation. These coefficients determine one disorder 
realization. In addition to a linear chain of qubits we also 
analyzed a case with qubits distributed on a square lattice 
which gave qualitatively similar results (see Sec. 6). 



4 Quantum computation with random errors 

The numerical results for fidelity decay induced by ran- 
dom errors in quantum gates are presented in Fig.lHl They 
clearly show that the decrease is exponential with time t 
and is given by the fit: 

f{t) = exp{-t/tr) ; tr = l/(0.095£^n^) ^ 47/(£^n^) . 

(22) 

As discussed in the introduction, the decay rate per gate is 
proportional to since on each step noise transfer such a 
probability from an ideal state to all other states (see [23 
[T?)ll 13111 41136113 7 ' ) . With a few percent accuracy the numer- 
ical constant in ((22j) is close to the lower bound discussed 
in inZI. 

While with random errors the fidelity drops by a sig- 
nificant amount in a purely exponential way the situation 
in the case of static imperfections is more complicated (see 
Fig-inj- In this case the initial exponential decrease is fol- 
lowed by a gaussian exponential one. As a result static 
imperfections give a faster decay of fidelity. The transi- 
tion between these two types of decay depends on the 
strength of imperfections e. Moreover, in scaled variables 
the weaker is s the stronger is the gaussian decrease of 
fidelity (see Fig[2l). We shall describe these phenomenon 
using the following RMT approach. 



5 Quantum computation with static 
imperfections: RMT approach 

Let us denote by U the unitary operator for the quan- 
tum map with static imperfections and by U the unitary 
operator for the ideal quantum map. We denote by Uj^ 
j = 1, . . . the elementary quantum gates which con- 
stitute the quantum map. According to the description of 
the quantum algorithm in section |S1 we write : 



U = Un-...-U2-Ui 



and 



ASH 



(23) 



(24) 



where 5H is the hermitian operator describing the static 
imperfections. In numerical simulations we have used the 
particular expression ^21^ for 5H but we mention that 
our approach does not rely on this expression and is much 
more general. We now introduce an effective perturbation 
operator for the full quantum map by: U = U e^^^""^^ . The 
operator ^i^eff is determined by: 



(25) 
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Fig. 3. Data points (d) represent the fidehty decay for the 
case of random noise errors in the gates of quantum algorithm. 
The fidehty / = |<'0(t)|'0noise(t)>|^ is shown as a function of 
the scahng variable x = 0.095 ts^n^ with Uq = 10, 12, 14, 16 
and 0.001 < £ < 0.1. The full line corresponds to the function 
/ = exp(— x). Data points (a), (b) and (c) represent the fidelity 
decay (/ = \<ip(t)\ipst&t{t)>\'^) as a function of the scaling 
variable x = te^ riqui for the case of static imperfections with 



10, Tig 



2 '^(1 



+ 4 being the number of elementary 



gates for one application of the quantum map and with (a) 
£ = 3 • 10-^ (b) s = 6 • 10-^ and (c) e = 5 • 10"^. The scaling 
variable x corresponds in all cases to Ft where F is the decay 
rate obtained from Fermi's golden rule. Therefore, the data 
curves (a), (b), (c) are tangent to exp(— x) for small x. The 
number of shown data points has been strongly reduced in 
order to increase the visibility of the different symbols, (same 
approach is used for other figures). Here and in Figs. 4-8, the 
initial state is as in Fig. 2, K ^ 1.7, T = 27r/A^. 



with 



5H{3) 



C/f ^ 5H Ui-...- Uj-i 



(26) 



We mention that the precise relation between ^i^eff and 
6H is not important for the following argumentation and 
we will need only one characteristic time scale tc defined 
by: 

2^ (27) 



= ^tr (5ifeff 



We furthermore assume that tr((5iJeff ) = (the case of 
tr((5iJeff) 7^ can be trivially transformed to the this case 
PH and 6H in Eq. has actually a vanishing trace). 

Following Ref. 31 , we express the fidelity in terms of a 
correlation function of the perturbation. For this we write 
the fidehty at time t as f{t) = \ A{t)\'^ with the amphtude: 



A{t) = ([/-* {Ue'^"'''') 



,i5H^{t(t-r) . 



and 



(28) 



(29) 



Here (• • •)q denotes the quantum expectation value. For 
a fixed initial state |'0o> this expectation value is given 
by* (' ' ')q — ^V^ol • • • |'0o>- However, in the following we 
average over all possible initial states that corresponds to: 
(• • •)q = -^tr(...). Since we are interested in the case 
where the fidelity is close to 1 we can expand (|28j) up to 
second order in ^i^eff (or equivalently in e) : 

A{t) « 1 + i ^ (<5J?eff (t))q - - ^ {5H,^{t))q 

r=0 r=0 
t-1 n-l 

"2 E E ('^^eff(ri)(5ffeff(r2))Q . (30) 

Ti=0 T2=0 

Now we introduce the correlation function C(r) by: 

C(ti - T2) = U {SHMri)SH^s{T2))Q (31) 

= te (<5ifeff(Tl - r2)<5J?eff(0))o 



and we note that {5Hes{T)) q = j^tr{SHes) 
bining Eqs. (jUI), (POj) and we obtain: 



fit) « 1 



t 2^, 

r-rE(* 

r = l 



■r)C(r). 



0. Com- 



(32) 



This provides the general expression relating the fidelity 
and the correlation function (pTT|) previously obtained in 
Ref. 1ST]. 

We now assume that the unitary quantum map U can 
be modeled by a random matrix drawn from Dyson's cir- 
cular orthogonal (/? = 1) or unitary {P = 2) ensemble 
39 ,40 . As we will see it is useful to express U in terms of 
its eigenvectors and eigenphases: 



\ 



(33) 



On 



The matrix V is either real orthogonal {/3 = 1) or com- 
plex unitary (/3 = 2). Inserting Eqs. (jSHI) in (|2B), we 
obtain: 



(34) 



In the following, we want to evaluate the average of C(r) 
with respect to U . We first evaluate the average with re- 
spect to the matrix elements of V which gives for 1 : 



{C{t))u = ( ^ 



1 ^/v 



(35) 



\j,k=i 



Here we have used that tr((5i^eff) = and we have re- 
placed tr(Ji^eff^) = N/tc according to Eq. ((77|) . We have 
furthermore neglected corrections of order l/A/"^ which 
arise from small correlations between matrix elements of 
V at different positions. We note that in the first 
term vanishes for = 2. For P = 1 this term arises from 
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additional contributions (in the V-average) because the 
elements of V are real for this case. The diagonal contri- 
butions with j = k in the second term of ((35j) provide 
the constant (which simplifies the first term). The 
average over the non-diagonal contributions with j ^ k 
can be expressed |39E0] in terms of a double integral over 
the two-point density for the eigenphases Oj. Since this 
two-point density is related to the two-point correlation 
function of the random matrix theory we obtain for r > 1 : 



where 



1 

TV 



ds Y2{s) e 



27ri f s 



(36) 



(37) 



is the "two-level form factor" defined as the Fourier trans- 
form of the two-point correlation function Y2{s). The form 
factor of the Wigner-Dyson ensembles is well known 38 , 
IS^lEHj . For the unitary and orthogonal symmetry class it 
reads (in the large N limit): 



P-l . b{T)-< Q if|f|>l, 



(38) 



ri-2|f| + |f|ln(2|f| + l) if |f| <1, 
■ = { -l + |f|ln(|^) if|f|>l. 

(39) 

Inserting the average correlation function (|36|) in (|32|) and 
replacing the discrete sum by an integral, we finally obtain 
the following scaling expression for the fidelity: 



with 



(inm).-fx(l^ 



^X(<5 



) = -2 / df{s-f)h2(f) 
Jo 



(40) 

(41) 

(42) 



where s = t/N. Using the random matrix expressions (|H8|) . 
we find for (3 = 2: 




if 5 > 1 . 



(43) 



^(-3s-24s2 + 17s3) 

+ 3^(l + 3s-4s3)ln(2s+l) 



|ln(3)(s 



-i(2-3s + s2) 
-^(l-3s + 4s3)ln(2s 
-j^(l + 3s-4s3)ln(2s 



1) 
1) 



if s < 1 , 



> if s > 1 . 



(44) 



Eqs. (|4()I44|) provide the key result of this section. From 
the practical point of view the contribution of Sx{s) in 
((^T|) is not very important, since (for P = 1): 



-s^ + ^s^ (if 5 <C1) , 



(45) 



Sx{s) « ln(3) - l) s + ^ ln(2s) - ^ ln(3) + ^ 

« -0.17604s + - ln(2s) - 0.49063 
6 

(if s > 1) . 



(46) 



Therefore for small s the linear term and for large s the 
quadratic term dominate the behavior of xi^) the ex- 
pression (|4H) . In the next section we compare this theo- 
retical random matrix result to the numerical data of the 
fidelity obtained for the quantum version of the tent map. 

Before doing so, we want to discuss three particular 
points. First, we have to evaluate the time scale tc that 
characterizes the effective strength of the perturbation. 
From Eqs. (|^^ and (|27|) . we obtain in lowest order in s: 

T = N^ tr(.5i7(j)Ji7(fc)) (47) 

with SH{j) given by (|26j) . This expression is similar in 
structure to Eqs. ((Sn|) or but with the important 

difference that here the "time" index corresponds to the 
number of elementary gates and not to the iteration num- 
ber of the quantum map. Since the elementary gates affect 
only one or two qubits ( "spins" ) the correlation decay be- 
tween SH{j) and 5H{k) will be quite weak and we can 
obtain a good estimate of ((Tf)) by: 



tr. 



2 r 2 

arig 5 UqS 



(48) 



where a < 1 is a numerical constant taking into account 
the exact correlation decay and the trace has been eval- 
uated using Eq. ((2U. The numerical results of the next 
Section indicate clearly that a ^ 1/5 such that the overall 
numerical factor is 1 and we have: 



1 



9 2 
2^^ 



11 



4 . 



(49) 



The second point to discuss concerns the fact that the 
phase space is mixed and not completely chaotic. As it can 
be seen from Fig.[2lthe chaotic region fills approximately a 
fraction of 0.65 of the full phase space. If the initial state is 
a gaussian wave packet placed in the chaotic region then 
its penetration inside the integrable islands induced by 
quantum tunneling will take exponentially long time scale: 
(X exp(A/'). Therefore with a good approximation we may 
say that in absence of imperfections the dynamics takes 
place only inside the chaotic component. Let us introduce 
now the Heisenberg time scale = 2"^^^ which is deter- 
mined by an average energy level-spacing for 2^-? quantum 
levels in the whole phase space. If the whole phase space 
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is chaotic then in the above RMT approach = • How- 
ever, for the case of Fig. [2l the chaotic component covers 
only a relative fraction <j = 0.65 of the whole phase space. 
Due to that in the previous RMT expressions we should 
put N ^ atu = 0.65 • to determine properly the num- 
ber of chaotic states. As a consequence, the expression 
((lOj) now reads (/3 = 1 for the tent map): 



crtu 



-X 



t 



t 



2 



tc tctu 



(50) 



Here we have neglected the contribution of Sx{s). 

It is interesting to note that there is a formal analogy 
between the fidelity decay given by ((5n|) and the decrease 
of the probability to stay near the origin which has been 
studied in mesoscopic and RMT systems (see e.g. 49,50, 
ICTl52ll5Hp . There, the time scale tc is replaced by the dif- 
fusive Thouless time scale txh and the second term with 
the Heisenberg time scale has negative sign. 

Finally, the third point concerns the fact that the above 
results are based on the assumption that the quantum evo- 
lution given by the exact quantum algorithm can be de- 
scribed by RMT. In particular, we assume that in the dy- 
namical evolution the ergodicity is established very rapidly 
after a few map iterations. This is correct for the choice 
T = 27t/N which corresponds to the dynamics in one clas- 
sical cell. We note that it is also possible to have many 
classical cells by the alternative choice T = 2ttL/N with 
L ^ 1 but fixed in the semiclassical limit N ^ oo. For 
K above the global chaos border, the classical dynam- 
ics is governed by a diffusive dynamics which covers all 
cells after the Thouless time scale txh ^ N'^ / D ~ / 
where D is the diffusion constant given by Eq. (0). In this 
case the theoretical treatment has to be modified since the 
matrix U will not be a member of the circular ensemble. 
However, choosing a static perturbation sufficiently com- 
plicated such that it can be modeled by a random matrix, 
one can show that the relations ((l0j) - (|42|) relating the fi- 
delity to the two-level form factor are still valid. The two- 
point energy level correlation function for diffusive metals 
has been calculated in the frame work of diagrammatic 
perturbation theory by Altshuler and Shklovskii 58 (see 
also the review [SIl)- The two-level form factor is now 
given by 62 (r) = &2,RM(r) + 62,diff.(^) where 62,rm denotes 
the random matrix expressions ((38j) or ((39j) and 62,diff.(^) 
is the correction due to diffusive dynamics which is given 
for a cubic sample by 



^2,diff.(^) 



2 



(51) 



Here d is the spatial dimension and g the dimensionless 
conductance with g ~ tn/^xh and txh being the diffusive 
Thouless time. In the limit f <^ g~^ (corresponding to: 
t <^ txh since f = t/N with N ^ tu) the sum can be 
approximated by an integral: 



y2,diff 



(t) 



2 fi-''/^ 



(52) 



Inserting this in we obtain the following diffusive 

correction to the scaling function ((^T|) : 



(3{2-d/2){3-d/2) {47rg)d/^ 



(53) 



We note that this contribution dominates the nonlinear 
RMT correction to x{s) in (gll) for d = 1, 2, 3 if = 2 and 
for (i = 3 if /3 = 1. The fidelity itself is slightly reduced by 
the diffusive correction according to : 



f{t) = /rm(^) exp 



4.3-d/2+d/2^ 



(54) 



where 5 is a positive numerical constant of order one. We 
mention that this interesting signature of the Altshuler- 
Shklovskii corrections for diffusive quantum systems in the 
fidelity decay is in principle accessible to efficient quantum 
computation. For our case with /3 = 1 and d = 1 this cor- 
rection is small. However for general quantum algorithms 
with diffusive behavior it may be important. The fact that 
this correction reduces the fidelity agrees with the obser- 
vation that the reduction of the volume of the chaotic 
component (a) also leads to faster fidelity decay accord- 
ing to Eq. (IHn|) . 



6 Quantum computation with static 
imperfections: numerical results 




0.0001 



1 10 100 1000 10000 

t 

Fig. 4. The fidelity decay for static imperfections [curves 
(a), (b), (c) with same values for riq and s as in Fig. 0] and 
random errors [curve (d) with riq = 10 and s = 1.59 • 10~^] as 
a function of t in a double logarithmic representation. The full 
lines correspond to power laws: — ln(/) ^ t and: — ln(/) ^ . 
The value of e for (d) is chosen such that the average reduction 
of fidelity for one elementary quantum gate is the same as for 
(c). The vertical dashed line provides the approximate position 
0.5crtH where the curves (a), (b), (c) change from linear to 
quadratic behavior. 
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We now consider the precise model of static imper- 
fections given by Eq. (|2T|) . We have numerically calcu- 
lated the fidelity f{t) for the tent map with K = 1.7 for 
Uq e {6, 8, 10, 12, 14, 16, 18} and 5 • lO"'^ < s < lO"^. 
For most cases we have determined the fidelity decay up 
to time scales t < tmax with f(tmax) = 0.5 (except for 
ng = 18 and the smallest values of e) since we are mostly 
interested in the regime (1 — /) <C 1 for which the an- 
alytical theory of the previous section is valid. We have 
also considered values £ > 10~^ but here the value tmax 
is typically so small that the number of available data 
points is not useful for the scaling analysis given below. 
In most cases we have concentrated on one particular re- 
alization of the random coefficients 5i and J^. But we also 
have made checks with up to 200 particular realizations. 
As initial state = 0)> we have chosen a coherent 
state \(p{po^Oo)> [see next section, Eqs. (|60|) . (|6T]) ] which 
is quite well localized around a classical point (po,^o) in 
phase space with a relative width ~ l/y/N - 2-^^/2 in 
both directions. 

First, we chose a state close to the hyperbolic fix point 
6> = 7r/2,p = 0, well inside the chaotic region of phase 
space. As can be seen in Fig. |21 after t = lb iterations the 
state fills up a big fraction of the chaotic region and after 
t 30 the state is practically ergodic. It covers then a 
fraction a ^ 0.65 of phase space. 




0.01 1 100 

Fig. 5. Scaling representation of the fidelity / for one particu- 
lar realization of static imperfections. The upper scaling curve 
shows: — ln(/)tc/tH as a function of t/tn with the two theoret- 
ical time scales tc — {e^nqn^g)~^ and tn — 2^"^. The full line in 
the upper curve corresponds to the theoretical random matrix 
result (|4()|) . (|41|) for (3 — 1. The lower scaling curve (shifted 
down by a factor 0.01) correspond to — ln(/)tc/tH versus t/tn 
with the times scales ic and tn obtained from the fit (us- 
ing appropriate weight- factors; see text) for each value of riq 
and £. Here the full line corresponds to the analytical scaling 
curve: y = x + (with y = — ln(/) tc/tn and x = t/tn). Data 
points are shown for s = 5 • 10~^ (and 6 < rig < 18); points 
for other values of s with 5 • 10~^ < ^ < 10~^ fall on the same 
scaling curve and are not shown. 



We have already seen in Fig. IHl of section El that the 
fidelity decay for static imperfections is faster than the 
exponential behavior for random errors. In order to an- 
alyze this in more detail we show in a double logarith- 
mic representation in Fig. ^ — ln(/(t)) as a function of 
t for the three cases already shown in Fig. [2l {uq = 10 
and e = 3 ' 10~^, 6 • 10~^, 5 • 10~^). For comparison, we 
also provide one case for random errors {uq = 10 and 
e = 1.59-10-'^). 

For the static imperfections, we can clearly identify 
a transition from linear to quadratic behavior at a time 
scale O.batu ^ 0.325tH corresponding to the theoretical 
expression ((SOI). However, the quadratic regime is best 
visible for the smallest values of s due to the restriction 
— ln(/) < ln(2). For the case of random errors there is no 
such transition and the linear behavior applies for all time 
scales. 

To analyze this transition in a quantitative way we 
determine for each value of Uq and s two time scales ic 
and in by the numerical fit 

y(t) = l + jL , y{t) = -Hf{t)). (55) 

tc tctu 

In order to prevent this fit to be artificially dominated by 
the large values of t (i.e. the quadratic regime) we mini- 
mize: 

d{ao, «i) = ^ ^(^) [y{t) - aot - ait'^] ^ (56) 
t 

with an appropriate weight factor w{t) ~ l/{ty'^{t)). The 
factor 1/y^ ensures that the vertical distance to be mini- 
mized is measured in the logarithmic representation for y. 
The other factor 1/t takes into account that the horizon- 
tal density of data points in the logarithmic representa- 
tion increases with t. The fit procedure (|56|) corresponds 
therefore to a fit in log-log representation such that also 
the small time scales (and values of y) are taken properly 
into account. 

According to the theoretical expression (|50|) one ex- 
pects that: 

ic = tc= o ^ o ^ tH = 0.5 cr tH ~ 0.325 tu (57) 

S^UqU^g 

with = 2^*2 and Ug = — ^ riq + 4. This theoretical 
prediction is verified in Figs. I5l8i 

In Fig. 13 we show two types of scaling curves for 
the fidelity. The first (upper) curve shows: — ln(/)tc/tH 
versus t/t^ with the time scales tc and given above. 
We observe that the numerical data coincide very well 
for Uq > 10 with the analytical random matrix result 
((l0j) . (|41|) for /3 = 1. The data for = 6, 8 show a moder- 
ate deviation for t > t^. We note that for this first scaling 
curve the dependence of the scaling parameters tc and 
on e and Uq is entirely determined by their theoretical ex- 
pressions. This is different for the second (lower) scaling 
curve where: — ln(/)tc/?H versus t/iu is shown. Here the 
scaling parameters ic and i^ have been obtained by the fit 
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for each value of Uq and e. Therefore all data coincide 
well with analytical scaling expression (|55j) . 

It is important to note that both scaling curves cover 
10 orders of magnitude and provide a strong confirma- 
tion of the crossover from linear to quadratic behavior 
predicted by the RMT approach. We mention as a side 
remark that we have also performed a similar scaling anal- 
ysis for the case of random errors. Here the scaling curve 
is purely linear in accordance with Fig. [2l However, this 
gives a stronger confirmation of the linear behavior than 
in Fig. ini since there the data for small e and large Uq 
corresponding to the regime (1 — /) <C 1 are quite badly 
visible in contrast to the scaling curve. 

In Figs. El and [3 the time scales ic and in obtained 
from the fit are shown versus tc and tn- We observe 
that the first theoretical expectation ic = tc is very well 
verified for the majority of data points. The small devia- 
tion for the remaining points appear for small Uq and the 
largest values of s where the fit procedure is less reliable. 
The second identity tu = 0.325 is in general also quite 
well verified. However, the deviations are slightly larger 
especially for larger values of £. Furthermore, for Uq > 14 
the regime t is numerically not accessible and the 

fit procedure amounts to extrapolate from data points 
t ^ tu or even t < tu for Uq = 16, 18. 

We also note that the data points for = 6, 8 (lower 
panel of Fig.EJ lie above the theoretical line in accordance 
with the first scaling curve in Fig. 

We have also determined the time scale at which 
f{tf) = 0.9. The theoretical expression suggests : 

^ _ 2tc ln(f ) _ 0.2107 -tc 

(58) 

For tc <C this implies tf ^ tc while for tc ^ 

we have tf ^ \/t^t^ = \/t~c2^^^'^. In Fig. [HI we show tf ob- 
tained from the numerical data for £ = 10"^, 5-10~^ and 
all values of Uq. The data points for > 10 coincide very 
well with while the points for = 6, 8 lie slightly 
above the theoretical line. For comparison, we also show 
the time scale t / obtained from the simplified exponential 
behavior f{t) = ex.p{—t/tc). Generally, tf is believed to 
decrease with increasing Uq and fixed s. However, this is 
not the case if tc > tn, i.e. for e < 2~^«/^/(n^y^). For 
very small values of e there is certain regime where t / first 
slightly increases with Uq and then decreases. 

In all presented numerical studies we considered the 
static couplings between qubits ordered on a line. To check 
that the results are not sensitive to this specific configura- 
tion we also considered the case when qubits are located 
on a square lattice as it was discussed in 19 . The ob- 
tained results (that we do not show here) confirms the 
RMT scaling (pjl . 

The case of the quantum evolution inside the inte- 
grable component of the tent map is analyzed in Figs. El- 
im Here, the initial state is located at ^ = 5.35 and p = 
which is in middle between the center fix point (6> = 37r/2, 
p = 0) and the boundary of the stable island {0 ~ 6.0, 
p = 0). We have determined for this case the time scales 




10^ 10^ lo'^ 

tc 

Fig. 6. The time scale tc obtained from the fit: — ln(/) = 

t/ic + tyiiciu) versus the theoretical expression tc — 

{e^nqn^g)~^ in a double logarithmic representation. The full 
line corresponds to ic = tc. The data points of top panel corre- 
spond to the same realization of static imperfections as in Fig. 
Elwith 6 < < 18 and 5- 10"'^ <s < 10"^. The bottom panel 
shows again the data points for 6 < riq < 18 and s = 5 • 10~^ 
(a) and e = 10"^ (b) for the same realization. The data points 
(c) are obtained from (t~^)~"^ where (• • •) denotes the average 
over 200 realizations of static imperfections for each value of 
Uq — 6, 8, 10, 12, 14 and e = 5 • 10~^ (statistical error bars are 
smaller than symbol size). 



ic and from the fit and performed the same scaling 
analysis in Fig. [HI for the fidelity decay as in in Fig. for 
the initial condition in the chaotic component. The scaling 
curves with the theoretical expressions (|F7)) for tc and 
give a significant deviation from the RMT result (upper 
group of curves in Fig. [2|). To understand the reason of 
this dispersion we also show the scaling curves with the 
fitted time scales ic and i^. This procedure gives a good 
scaling of numerical data (lower group of curves in Fig. El). 
Obviously, the fit (|55j) still works quite well as such but the 
obtained fit parameters are eventually different from the 
initial condition in the chaotic component and the RMT. 
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Fig. 7. The time scale tn obtained from the fit: — ln(/) = 
versus the Heisenberg time tn = 2^*2 in a 
double logarithmic representation. The full line corresponds 
to in = O.batu- The data points in the top panel correspond 
to the same realization of static imperfections as in Fig. m 
with 6 < < 18 and 5 • 10"'^ < £ < 10~^. The bottom 
panel shows again the data points for 6 < rig < 18 and 
£ = 5 • 10~^ [data points (a)] and e — 10~^ [data points 
(b)] for the same realization. The data points (c) are ob- 
tained from {ic^){{tciB)~^)~^ where (• • •) denotes the average 
over 200 realizations of static imperfections for each value of 
riq — 6, 8, 10, 12, 14 and £ = 5 • 10~^ (statistical error bars are 
smaller than symbol size). 



The dependence of tc on the theoretical value of tc is 
presented in Fig. It shows that the theoretical expres- 
sion works with a good accuracy in the interval of 6 orders 
of magnitude. This is not really a surprise since according 
to ^ t-^ measures the overall strength of the pertur- 
bation. However, for shown in Fig. ^2 the situation is 
much more complicated. The variation of vs. = 2^*2 
shows unusual steps and it is unclear what is the real 
dependence in the limit of large Uq. Further studies are 
required for complete understanding of the static imper- 
fection effects in this regime. This fact compromises the 



Fig. 8. The time scale tf determined by /(t/) = 0.9 (in a log- 
arithmic representation) as a function of the number of qubits 
Tiq. The data points correspond to the numerical simulation 
for the same realization of static imperfections as in Fig. m 
with £ = 5 • 10"'^ (a) and e = 10"^ (b). The two full lines 
correspond to t/ given by assuming the theoretical ex- 
pression Eq. (|50f) for the fidelity. The two dashed curves corre- 
spond to t/ = tc In(^) for the simplified exponential behavior 
/(t)=exp(-t/te). 



possibility to determine the asymptotic dependence of t/ 
on £ and Uq for the case of integrable or quasi-integrable 
dynamics. In addition, the data of Fig.|Hlshow that it is not 
easy to determine the asymptotic behavior oft/ in absence 
of clear scaling laws. Due to these two remarks we think 
that the scaling dependence for tf time scale, proposed in 
[T51l55j for the case of static imperfections, represents in 
fact only an intermediate behavior and cannot be extrap- 
olated to the limit of large Uq. Indeed, the quantum evo- 
lutions studied in 13,36 correspond to quasi-integrable 
regimes and additional tests are required to check the va- 
lidity of the RMT scaling (|50j) for the quantum algorithms 
studied there. 

Finally, it is interesting to compare directly the fi- 
delity decay induced by static imperfections for the quan- 
tum evolution in chaotic and integrable components (see 
Fig [T2j) . The numerical data show that fit) decreases faster 
in the case of integrable evolution. As it was discussed in 
pTT) the presence of chaos reduces the fidelity decay rate. 
This is in the agreement with the results of Fig. [3 and 
Fig. ^3 according to which tn is much smaller for the in- 
tegrable regime as for the regime of quantum chaos. How- 
ever, the possibility of using this fact to improve the accu- 
racy of quantum computations remains an open question. 

The numerical data presented in this section definitely 
confirm the RMT universal law for fidelity decay for the 
case when the evolution takes place in the regime of quan- 
tum chaos. This means that this law works for quantum 
algorithms simulating a complex dynamics. The situation 
for the evolution in the integrable component is more com- 
plicated. The data show that the theoretical expression for 
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0.01 1 100 

Fig. 9. ScaUng curve for fidehty as in the case of Fig. but 
for the quantum evolution inside integrable component. Here, 
the initial state is a minimal coherent wave packet taken inside 
the regular part of phase space at ^ = 5.35 and p = 0; the same 
realization of static imperfections as in Fig. El is used. 
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Fig. 11. The time scale iu as in Fig. versus tu = 2""^ 
data are obtained from Fig.^jfor the case of regular dynamics. 
Shown are data points for e = 5 • 10~^ (a) and e = 10~^ (b) 
with 6 < rig < 18 in a double logarithmic representation. The 
full line corresponds to tn = 0.5crtH- 
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Fig. 10. The time scale tc as in Fig.Elversus tc = (e^rign^)""^, 
data are obtained from Fig.^lfor the case of regular dynamics. 
Shown are data points for e = 5 • 10~^ (a) and e — 10~^ (b) 
with 6 < rig < 18 in a double logarithmic representation. The 
full line corresponds to tc — te- 
le is still valid but the dependence of the time scale on 
Tiq requires further investigations. 

It is interesting to note that the relation (|50j) should 
also work for the problem of Loschmidt echo in systems 
with quantum chaos j27ll28llH()pHl ' 32 . In this case for small 
perturbations tc is still given by Eq. (^7)) or, that is equiva- 
lent, the inverse decrease rate is given by the Fermi golden 
rule [28|. Then the scale is determined by the inverse 
density of states or for quantum maps by the number of 
states via relation = A/". As a result for small per- 
turbations the decay of Loschmidt echo for such quan- 
tum dynamics is still given by the universal decay relation 
Eq. JSII). 




t 

Fig. 12. fidelity decay for initial minimal coherent state in 
chaotic (full curve) and integrable (dashed curve) component 
for rig = 16, £ = 5 • 10~^. 

7 Husimi function 

Here we discuss how an arbitrary quantum state \\\)> can 
be represented in the classical phase space in the process of 
quantum computation. For this it is convenient to use the 
coarse-grained Wigner function (or the Husimi function) 

9h{.vM = \<^^VM\^>? (59) 
where the smoothing is done with the coherent state 

|(^(po,^o)> = A^e-^^-^o)'/^'^'-^^^^ \v> . (60) 
p 

Here, A is the normalization constant and a is the width 
of the coherent state in the ^-representation. The coherent 
state corresponds to a gaussian wave packet that is local- 
ized in the classical phase space around a point (^o, Po) 
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with widths Ap = a and AO = l/(2a). We choose a = 
such that the widths relative to the size of the 
phase space are comparable: 

Ap _ 1 _ 0.2887 AO _ ^ 0.2757 

N ~ VT27V ' 27r ~ 27rVN ^ VN 

(61) 

The naive evaluation of the Husimi function with- 
out any optimization requires 0{NNpNo) operations (on 
a classical computer) where Np and Nq are the numbers of 
values for po and for which is evaluated. In view of 
Eq. (|61|) it is sufficient to choose Np = Nq = ^/N resulting 
in 0{N'^) operations which is very expensive as compared 
to 0{N\og{N)'^) operations needed by the simulation of 
the quantum map on a classical computer as described in 
Sec. 01 

Fortunately, the evaluation of the Husimi function can 
be done in a more efficient way. To motivate and explain 
this let us first consider a modified Husimi function defined 
by: 

p^p{po,0o) = \<^^''Kpo,0o)\4'>\' (62) 
with the modified coherent state: 

\^ip)(p„eo)> = ^ E e"''"^b>- (63) 

^ P=Po 

Here we assume for the sake of simplicity that the number 
of qubits Uq is even such that a/TV = 2^^/^ is integer. We 
furthermore require that po is an integer multiple of a/TV 
and Oo 2tt 1/\/N with / = 0, . . . , VTV - 1. 

Comparing with ((^ . we see that the gaussian 
pre-factor has been replaced by a box-function of width 
a/TV. This provides a very good localization for the mo- 
mentum representation but implies that in angle repre- 
sentation the amplitude around Oq decreases only as a 
power law according to the Fourier transform of the box- 
function: sin(x)/x with x = ^/N{0 — 0o)/2. However, the 
modified coherent state ((63j) still provides a quite well lo- 
calized state around the point (^o, Po) - Its main advantage 
is related to the fact that it can be put in the form: 

\v>^''Hpo,0o)> = Uqft\Po + 1> (64) 



first half of the binary digits the information for and in 
its second half the information for po. More explicitly, if 
p = Yl^lo^ aj2^ , we have : 

riq-l ng/2-1 

Po= ^^-2' ' ^0 = ^ E ^^-2' • (66) 

We note that it is also possible to introduce another 
type of modified Husimi function (and modified coherent 
state) by exchanging the roles of Oq and po : 

pfiPo, Oo) = \<P I C/q^t f^QFT I i'>\^ ■ (67) 

where p = N0o/{27r) -\-po/\/N. As in Eq. (|?)^ . we require 
that Po is an integer multiple of and Oq = 2tiI/^/N 
with / = 0, . . . , \fN — 1. The operator /7qft corresponds 
to quantum Fourier transform for all qubits and trans- 
forms a state from p- to 6>-representation. [/qft corre- 
sponds as above to quantum Fourier transform for the first 
half of the qubits. We mention that the coherent states as- 
sociated to ((S7)) have a power-law localization amplitude 
for p and a box-function localization amplitude for 0. 

We have seen that both types of modified Husimi 
functions can be evaluated on a classical computer with 
0{N\og{NY) operations. Based on the idea of the QFT 
which is closely related to the FFT if simulated on a clas- 
sical computer, we have also implemented an efficient clas- 
sical algorithm for the original Husimi function ((5^ with 
the gaussian coherent states. For each value ofpo, we only 
consider the restricted sum such that \p — po\ < 4a/]V and 
evaluate the matrix elements <(p{po^ Oq) \ iIj> for all values 
of ^0 simultaneously using FFT. This provides also an al- 
gorithm with complexity 0{N \og{N)'^) but with a consid- 
erably larger numerical pre-factor. However, this method 
does not allow for a "pure" quantum computation as it is 
possible for the two types of modified Husimi functions. 

In order to compare the different Husimi functions, we 
consider a test-state defined by a circular superposition of 
coherent states as: 

|V^circle> =A \ip{po, Oo)> • (68) 



where Oq = 2tiI/^/N and /7qft corresponds to the quan- Here i is a normalization constant and the sum runs over 

t um Fourier transform operator (see (HSI)) for the first half ^ discrete set of points (i9o,Po) on a circle with center 

of the qubits (ao, Q^n,/2-i). Eq. ^ implies that (tt, A^/2) and relative diameter 0.7 (as compared to the 

size of the phase space). 

Pi? (Po,^o) = \<p\Uqyt I V^>P (65) In Fig. [TSl we show the density plots for the three 

types of Husimi functions for this test-state with Uq = 

withp = po-^l = po^VN0o/{27r). Here the state L^qft | 10, 12, 14. In all cases the circle picture of the density is 

can be evaluated efficiently on a quantum computer us- quite well reproduced and one clearly sees that the circle 

ing + 1) elementary quantum gates according to has a finite width according to the widths of the coher- 

(Uni) (with Uq replaced by nq/2). Emulating the quan- ent states due to the quantum uncertainty principle (see 

tum computer on a classical computer this still costs only Eq. ^^). For the modified Husimi function (column 

0{N\og{N)'^) elementary operations. The matrix elements (b)), one clearly observes the effect of the power-law de- 

<p\ Uqft I '0> of this state with the momentum eigen- crease for the ^-amplitude leading to a smearing out of 

states \p> provide directly via Eq. the modified Husimi maxima in the 6>-direction. The same holds for the second 

function. Here the value of p = 0, ... A^ — 1 contains in its modified Husimi function (|?)7)) (column (c)) concerning 
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Fig. 13. Density plots of the three different types of Husimi 
functions for the circle-state with Uq = 10, 12, 14 

(rows from top to bottom). The horizontal axis corresponds 
to ^0 ^ [0, 27r[ and the vertical axis to po G [0, A/'[. The density 
is minimal for blue/black and maximal for red/white. Column 
(a) corresponds to the Husimi function (|59|) with gaussian am- 
plitude, column (b) to the momentum modified Husimi func- 
tion and column (c) to the phase modified Husimi function 



the j9-direction. This effect is strongest for small values of 
riq and becomes smaller with increasing Uq. The effect of 
smearing out is not visible for the original Husimi function 
(|59j) (column (a)) with gaussian amplitudes for 6 and p. 

In order, to study the evolution of the Husimi func- 
tion for chaotic and regular regimes in the quantum tent 
map, we choose as initial condition the circle-state (|68|) . 
The semi-classical density of this state intersects quite well 
with both regular and chaotic parts of the mixed phase 
space (see Figs. I1I2| . 

In Fig. ^2 we show the density-plots of the Husimi 
functions (defined by Eq. of the state obtained from 
the circle-state after 100 iterations of the quantum map for 
the cases Uq = 10, 12, 14 (column (a)). We also show the 
state that is obtained by applying further 100 iterations 
of the inverse quantum map which should theoretically 
provide the original circle-state (column (c)). In Fig. IT^ 
we also show a density-plot for the classical map ^ (with 
p to be taken modulo 27r/T). Here we have determined 
the classical trajectories of lOON random initial points on 
the circle. Then the density-plot has been calculated from 
a histogram with a finite box-size corresponding to the 
finite resolution of the quantum case with Uq = 14. 



Fig. 14. Density plot of the Husimi function (|59|) from the 
the circle-state (j68|) after 100 iterations with the quantum map 
with K = kT = 1.7 and T = 27r/Ar = 27r/2^'? (columns (a) 
and (b)) and after 100 iterations in the future and 100 itera- 
tions with the inverse map in the past (columns (c) and (d)). 
The first three rows correspond to Uq — 10, 12, 14 (from top 
to bottom) and the last row corresponds a classical density 
plot obtained from a histogram-sampling with a box-size cor- 
responding to the resolution for Uq — 14 and with an average 
number of 100 classical trajectories per box. Columns (a) and 
(c) correspond to the exact quantum or classical maps only 
limited by the relative machine precision 10~^^. Columns (b) 
and (d) correspond to the quantum computation with random 
errors in quantum gates [e — 0.01) or perturbed the classical 
map perturbed by noise with e — 0.01 (see text). 



One clearly sees in column (a) that the chaotic part of 
the phase space is filled up ergodically while the piece of 
the circle intersecting the regular part of the phase space 
remains a connected line. Actually, this line rotates with a 
constant angular velocity around the center fix-point due 
to the local linear behavior of V {6) close to the fix-point. 

Concerning the states obtained after the back-iteration 
in time, which are shown in column (c), one observes that 
the inverse quantum map reproduces exactly the initial 
state while for the classical map only the pieces of the cir- 
cle belonging to the regular part of the phase space are 
reproduced. This is due to the finite machine precision 
(of 10~^^) together with the exponential instability in the 
chaotic part of the phase space. We have verified that for 
only 25 iterations, the circle is well reproduced in every 
part of the phase space. At 50 iterations the classical com- 
puter round-off errors already have significant effects but 
are not sufficient to create a uniform distribution in the 
chaotic region as it is the case with 100 iterations shown in 
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Fig. last row of colunin (c). This effect is corapletely 
absent in the quantum simulation. The information for 
the phase space distribution is encoded in the quantum 
state in such a way that it is not sensible to the round-off 
errors of the classical computer simulating the quantum 
algorithm for the tent map. 

To investigate this point in more detail, we have also 
performed a quantum simulation where all quantum gates 
are perturbed by random errors (see Sec. 3). The effects 
of this noisy perturbation can be seen in Fig. in the 
columns (b) (100 forward iterations of the circle-state) and 
(d) (100 forward and 100 backward iterations) where we 
have chosen e = 0.01. Concerning the quantum map, the 
noise reduces some-how the general quality of the pictures 
but it does not distinguish between chaotic and regular re- 
gions of the phase space. In particular in column (d), the 
circular density is quite well reproduced with some ad- 
ditional overall noise. Concerning the classical map, the 
circle-pieces in the regular region still remain closed lines 
but they acquire a finite width which increases in a diffu- 
sive way with the number of iterations. The circle-pieces 
in the chaotic region become very quickly mixed. Further- 
more, it is not possible to reproduce the initial circle in 
the chaotic region due to the exponential instability (last 
row of column (d)). We have verified that this effect is al- 
ready true for only 15 forward and 15 backward iterations 
if 6 = 0.01 (for the classical map the noise is introduced in 
the equation for momentum with an amplitude e = 0.01, 
Fig. [Til bottom d). 

We have also studied the effects of static imperfec- 
tions on the Husimi function evolution in the tent map. 
In Fig. El we show the results of static imperfections with 
e = 10~^ and Uq = 10, 12, 14. The initial state is again 
the circle-state and column (a) corresponds to the state 
after 100 forward iterations and column (b) to the state 
after 100 forward and 100 backward iterations. The effect 
is quite similar to the quantum computation with random 
errors (columns (b) and (d) of Fig. ITU) . The general quality 
of the pictures is reduced and there is no distinction be- 
tween regular and chaotic part of the phase space. Again, 
in column (b) the circular density is quite well reproduced 
with some additional overall noise. We should note that 
the static imperfections of strength e = 10~^ give pertur- 
bations in the Husimi function which are comparable with 
those in the case of random errors at e = 0.01. This shows 
that the static imperfections perturb the quantum com- 
putations in a stronger way comparing to random errors. 

Finally, we show in Fig. the modified Husimi func- 
tions ((62j) after 100 iterations applied to the initial circle- 
state again for the three cases Uq = 10, 12, 14. Column (a) 
shows the exact simulation, (b) the case of random errors 
{s = 0.01) and (c) the quantum map with static imper- 
fections {e = 10~^). We note that the smear ing-out effect 
discussed at the beginning of this Section (see Fig. 
is well visible for the case of the exact simulation, while 
it is not visible at all for the cases with random errors 
or static imperfections. Therefore, the utilization of the 
modified Husimi function seems to be quite well justified 
in these cases. 



a b 




Fig. 15. Density plot of the Husimi function (|59|) from the 
the circle-state (|68|) after 100 iterations with the quantum map 
(same values as in Fig. ll4f) perturbed by static errors (see text) 
with e = 10~^ (column (a)) and after 100 iterations in the 
future and 100 iterations with the inverse perturbed map in 
the past (column (b)). The different rows correspond to Uq — 
10, 12, 14 (from top to bottom). 

8 Conclusion 

The results obtained in this paper give a universal de- 
scription of fidelity decay in quantum algorithms simulat- 
ing complex dynamics on a realistic quantum computer 
with static imperfections. This decay is given by Eq. EOI 
which determines the time scale of reliable quantum 
computation with fidelity / > 0.9. According to Eg. IKHI 

if ^ tc/10 = l/(10£Xn^) ; Ng ^ l/(10£Xn^) (69) 

for > tc so that e > £ch = 2~"''?/^/(n^y^). Here, 
Ng = tfUg is the total number of gates which can be 
performed with fidelity / > 0.9. In this regime the static 
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Fig. 16. Density plot of the modified Husimi function (|62|) 
from the the circle-state (|68|) after 100 iterations with the quan- 
tum map (same values as in Fig. I^)- Column (a) corresponds 
to the exact quantum map, column (b) to the map with ran- 
dom errors in quantum gates with e = 0.01 and column (c) 
to the quantum map simulated with static imperfections with 
e = 10"^ The different rows correspond to Uq — 10, 12, 14 
(from top to bottom). 

errors act in a way similar to random noise errors even 
if their effect is stronger due to coherent accumulation of 
static errors inside a certain interval of the algorithm (one 
map iteration for the tent map). Indeed, for random errors 
in quantum gates the relation ((221) gives 

if ^ t./lO ^ h/{e^ng) ; Ng ^ b/e^. (70) 

We note that (fZOj) is in agreement with the result obtained 
for random errors in a very different quantum algorithm 
ISni and hence it is generic. Even if the dependence of Ng 
on £ in Eqs. ((69j) . (fZOj) is the same, the dependence on 
Uq is rather different. This difference should play an im- 
portant role for the quantum error correction codes which 
allow to perform the fault-tolerant quantum computation 
for the random error rate Pr ^ < 10~^ (see e.g. 
I!ni54ll55ll5f)j ). The fact that for random errors Ng is in- 
dependent of Uq while for static imperfections Ng drops 
strongly with Uq should significantly decrease the thresh- 
old for fault-tolerant quantum computation in presence of 
static imperfections. 

For tc < tu or e < Sch = ^""''^/^/(npy^) the time 
scale is given by the relation 

tf ^ 0.2Vt^ ^ ^^'^l^liheug^q) ; Ng ^ 2"^/V(5^V^)- 

(71) 



In this regime the effect of static imperfections is abso- 
lutely different from random noise errors. This regime may 
be dominant for up to 10 - 15 qubits. However, in the limit 
of large Uq ^ 10 it appears only in the limit of very small 
static imperfections and should not be very important for 
quantum computers with few tens of qubits. The transi- 
tion from the regime ((69j) to regime ^7l\ takes place for 

s>eeh = 2-"''/V(ng^A^) • (72) 

From the physical point of view this border can be in- 
terpreted as the quantum chaos border above which the 
static imperfections start to mix the energy levels of ideal 
quantum algorithm. The fact that this border drops expo- 
nentially with the number of qubits riq has been discussed 
in for a quantum algorithm for complex dynamics. 
Above Sch the effect of static imperfections becomes some- 
what similar to random errors. 

The results (|69|) and ((7T|) for the time scales of reliable 
quantum computation are based on the RMT approach 
and are universal for algorithms which simulate a com- 
plex dynamics, e.g. an evolution in the regime of quan- 
tum chaos. However, it is important to keep in mind that 
there are other types of algorithms where the evolution is 
rather regular, e.g. the Grover algorithm or integrable dy- 
namics. In such cases the asymptotic dependence of on 
Uq should be studied in more detail. It is not excluded that 
in such cases tu grows with Uq very slowly (see Fig. ITT]) 
or even may be independent of riq. In such situations the 
static imperfections will generate a very significant reduc- 
tion of the time scale of reliable quantum computation. 
In a sense our RMT result ((50j) gives the weakest form 
of fidelity decay in a realistic quantum computer with Uq 
qubits since the reduction of the chaotic component a ac- 
celerates this decay. 

The universal regime for fidelity decay in quantum 
computations established in this paper can also find other 
applications. For example it can appear in the decay of 
spin echo in interacting spin systems. 
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Upon completion of this manuscript a recent preprint 
of T. Gorin, T. Prosen, and T. H. Seligman 59j came to 
our attention where the relation between fidelity decay 
and two-level form factor has been established for an ab- 
stract Hamiltonian model with continuous time evolution 
and a perturbation given by an invariant random matrix. 
However, in this work only the regime of a small pertur- 
bation strength A <C 1 (A ~ ytii/tc) is studied which 
translates to tu <C tc with the dominant gaussian decay. 
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